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Abstract 

We investigate the effect of the phase difference of apphed fields on the dynamics 
of mutually coupled Josephson junction. The system desynchronizes for any value 
of applied phase difference and the dynamics even changes from chaotic to periodic 
motion for some values of applied phase difference. We report that by keeping 
the value of phase difference as vr, the system continues to be in periodic motion 
for a wide range of system parameter values which might be of great practical 
applications. 
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1 Introduction 

The study of the dynamics of Josephson junction (JJ) is of fundamental and 
experimental interest. The interaction of Josephson junctions with external 
fields have played an important role in the development of Josephson physics 
and its chaotic dynamics [1,2,3,4]. The existence of chaos in rf-biased Joseph- 
son junction has been verified through theory, numerical simulation and ex- 
periments [5]. The rf-biased junction is of practical importance because of its 
application as voltage standards, a situation where chaotic behavior is least 
desired [6]. Control of chaos is an active area of research [7] because of the 
many undesirable effects chaos brings in mechanical systems and other devices. 
It was shown to be possible to control chaos both theoretically and experimen- 
tally using different methods such as giving a feed back [8], application of a 
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weak periodic force [9] etc. The problem of controlling spatiotemporal chaotic 
pattern induced by an applied rf signal in a Joscphson junction has earlier 
been discussed [10]. By controlling chaos in rf-biased Josephson junctions it 
was shown that even in the presence of thermal noise they can be used as volt- 
age standards [11]. Suppression of temporal and spatio-temporal chaos allows 
complex systems to be operated in highly nonlinear regimes. This is a desir- 
able feature in many physical systems. By applying a small timc-dcpcndcnt 
modulation to a parameter, a chaotic system can be stabilzed. However in 
practical applications this method requires that the characteristic times of the 
system is not too short compared with the times of the feed back system. In 
the case of JJ oscillators the characteristic times of the dynamics response arc 
of the orders of few picoseconds which is too short for any electronic feedback 
control system. 

Since it was shown that chaotic systems could be synchronized by linking them 
with common signal [12] many works have been done in this direction because 
of its application in secure communication [13]. Synchronization in Joseph- 
son junction has been an interesting area of research [14,15,16]. The role of 
phase difference of applied sinusoidal fields in desynchronizing and suppress- 
ing chaos in duffing oscillators has been studied [17,18]. In the present work, 
we consider the effect of phase difference of the applied rf-ficlds on mutually 
coupled Josephson junctions. We discuss the equation of coupled JJ in section 
II and arrive at the dimensionless first order form of the equation of motion 
of the system. In section III the parameter range in which the system is syn- 
chronized is found and the effect of phase difference of the applied rf fields on 
synchronization is found. We analyze the dynamics of the Josephson junction 
after applying the phase difference. Then by fixing a particular phase value 
where the system becomes periodic we discuss the effect of other parameters 
on the junction. In section IV the results are discussed and the applications 
are mentioned. 



2 THEORY 



Josephson junction can be represented by a resistively and capacitively shunted 
junction (RCSJ) model and the dynamics of the system can be explored by 
writing the equation of motion [19]. The equation of a single Josephson junc- 
tion by this model can be written by solving Kirchoff's law as 

hc h d(j) I . , 

2^ + 2^ + ^^'^"'^ = ^'^^ + ^° 

where is the phase difference of the wave function across the junction, 
iQCOsiujt) is the driving rf - field and idc is the dc bias. The junction is char- 
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acterized by a critical current ic, capacitance C and normal resistance R. The 
coupled J J considered here consists of a pair of such junctions wired in parallel 
with a linking resistor Rs [15]. Such a system can be schematically given as in 
Fig(l) and the dynamical equations can be written as 

Rs 

1 1 ^ , 




Fig. 1. Schematic representation of a coupled Josephson junction connected in par- 
allel with a linking resistor Rs. 1 and 2 represents the applied fields. 
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where is the current flowing through the coupling resistor and is given as 



h 

2eR, 



d(j)i d(f)2 



dt' 



dt' 



(4) 



In order to express Eqs.(2) and (3) in dimensionless, normalized form we 
introduce the junction plasma frequencies cuji and U!j2 which are given by 
ujji — (2eici/^Ci)2 and Uj2 = {2eic2/hC2)^ and also normalized time scale as 
t — cujit'. The dimensionless damping parameter is defined as 



h 



Ri V 2eiciCi 



The dc bias current i'^^ and the rf amplitude i'^ are normalized to the critical 
current id- We also normalize the actual frequency cu to — cu/cuji and 
as = (Ri/Rg) p. For identical Josephson junctions, we can write Eqs.(2) and 
(3) as 



01 + + sin 01 = idc + io cos{Qt) — ag 0i - 02 

02 + /502 + sin 02 = idc + io cos(Qi + 9) — ag 02 - 0i 



(5) 
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Fig. 2. Maxima of the normalized voltage plotted against coupling strength for 
e = 0, P = 0.15, io = 0.7 and u = 0.6. 
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Fig. 3. Maxima of the difference in voltage against coupling strength a^. ^ = 0, 
P = 0.15, io = 0.7 and u = 0.6. 

It can be seen that the coupling arises as a natural consequence of the exchange 
of current through the resistor Rg and it depends on the differential voltage 
(■01 — 1^2)- For Josephson junction devices, phase derivatives are of central 



importance because they are proportional to junction voltages. In order to 
study the system numerically we write the above equations in the first order 
differential form as 



Eq. (6) is studied using fourth order Runge-Kutta method and the maxima of 
normalized voltage values are plotted to study the dynamics. It was observed 
that the system exhibits chaotic behavior for j3 = 0.15, io = 0.7, idc = 0.3 and 
Q = 0.6 as seen in Fig. 2. From Fig. 3 we can see that the difference in voltage 
becomes smaller as the coupling strength is increased. 



3 The effect of Phase Difference 

The effect of phase difference of the applied sinusoidal driving field on coupled 
periodically driven duffing oscillator has been studied [18] and it was observed 
that the system desynchronizes by the application of a phase difference. To 
study how the phase difference of the applied fields acts on JJ system we write 
Scj) — (f>i — 02 and Sjjj — ipi — ip2 and from Eq. (6) we get 

S(i) — S^, (7) 
Sip — —fiS^ — sin 01 + sin 02 — '^OisS^ + 2io sin + sin ^-^ . 

When ^ = the term 2iosin [fit + 6/2) sin {6/2) vanishes. We choose the val- 
ues of as such that the difference in voltage is negligible, ie., ipi ~ ip2- Now 
both and go to zero. Prom Fig. 3 we select our value of ag as 0.45 which 
almost satisfies this condition. From Eq. 7 we can see that even for small 
values of phase differences ^ and the system desynchronizes. Figs. 4(a) 
and 4(c) show that the system is synchronized and Figs. 4(b) and 4(d) show 
that the system is desynchronized by an applied phase difference of ^ = O.Itt. 
The level of mismatch of chaos synchronization can be given quantitatively by 
taking the similarity function S{t) as a time averaged difference between the 
variables ipi and ipi taken with time shift r [20] 



01=V'l, 

02=^2, 
-02 = -Plp2 



sin 01 + ioc + io cos(Oi) - ag [ipi - 1P2] , 



sin 02 + idc + io cos[(Ot) + 6] - ag [ip2 - ^1] ■ 



(6) 




mt + r) - Mt)?) 



(8) 
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Wc plot S{t) against r for different values of phase difference 6 as shown in Fig. 
5. It is observed that for = 0, the system is in complete synchronization. 
For a finite value of phase difference we observe that a minimum of S{tq) 
appears which indicates the existence of a certain phase difference between 
the interacting systems. S{to) is finite in these cases which means that in this 
regime the amplitudes are uncorrelated. 
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Fig. 4. a and c show the system is synchronized for ^ = 0, /? = 0.15, iq = 0.7, 
idc = 0.3, as = 0.45 and u = 0.6. b and d show the system desynchronized for a 
phase difference of 9 = O.Itt 

Now we vary the values of 9 from to 27r and plot the maxima of the difference 
in voltage against the phase. It can be observed from Fig. 6 that the system 
passes from chaotic to periodic motion for some values of 6. In the synchroniza- 
tion manifold, ie., when 0i = 02 and ipi = ip2 we can write P,/, = [0i + ^2] /2 
and = [ipi + 1P2] /2 and from Eq. (6) we get 



P4> = P^, (9) 
= -(3P^ - sin(P<^) + idc + io cos j cos + ^ j . 

Eq. (9) is equivalent to Eq. (1) and we can see that is replaced by io cos {0/2) 
and the phase of the driving field leads by 6/2 as a result of coupling. This 
might correspond to a parameter space where the system is periodic which 
explains the change in system dynamics. 
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Fig. 5. Similarity function S{t) versus r for different values of phase difference 6. 
Curve 1 is with phase difference ^ = 0, 2 for ^ = O.Itt and 3 for ^ = O.Stt. 

The variation of the voltage of one junction with the appUed phase difference 
is shown in Fig. 7. From the inset of Fig. 6 it can be seen that the system 
exhibits periodic window in the region where 9 — O.SAtt — OAtt. However in 
this region even a shght change in system parameter values would bring the 
system back to chaotic regime. For a phase difference of ^ = 0.957r — I.Stt 
the system exhibits periodic motion and the corresponding periodic difference 
in voltage and voltage plotted against time is shown in Figs. 8(c) and 8(d). 
Figs. 8(a) and 8(b) shows difference in voltage and voltage against time for a 
phase difference of ^ = 0. 

Fixing the phase difference between the driving fields as vr, the change in 
the response of the system to other parameter variations are studied. From 
Fig. 2 we can see that the system is in chaotic motion even for large values 
of couphng strength. However by the application of a phase difference of tt 
we observe that the system exhibits periodic motion from a coupling strength 
0.32 onwards as can be seen from Fig. 9. 

Now we fix the value of a as 0.45 and the amplitude of the driving rf field 
is changed from to 1. Without an applied phase difference we can see from 
Fig. 10 that the system exhibits chaotic motion from a value of 0.43 onwards 
with some periodic windows in between. However on the application of a phase 
difference the system stays in periodic state for a wide range of amplitude 
values which were chaotic earlier. Fig. 11 shows the response of the system 
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Fig. 6. Maxima of the difference in voltage is plotted against phase difference applied, 
e = 0-2Tr. as = 0.45, idc = 0.3, (3 = 0.15, ^o = 0.7 and u = 0.6. 



when the amplitude of the driving field is changed from to 1 with an applied 
phase difference of tt. 

The effect of phase on the applied dc bias on the system is also studied. For 
this all other parameter values were fixed and idc value is changed from to 
0.4. Here also we observe that the system continues to be in periodic motion 
for a large range of i^c values. The comparison can be obtained from Fig. 12 
and Fig. 13. Thus we show that the system exhibits periodic motion for a 
wide range of parameter values for an applied phase difTcrcncc between the 
driving fields. This might be of great practical importance when we consider 
Josephson junction devices like voltage standards. 

An important point to be noted is that the parameter values at which we 
apply phase difference is to be chosen carefully. If the values we choose is in a 
region where the difference in voltage (■^i — 1P2) is large, then by just applying 
a phase difference we may not be able to control chaos. 
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Fig. 7. Voltage corresponding to a single junction is plotted for phase values 
6* = - 27r with = 0.45, = 0.3, , (3 = 0.15, iq = 0.7 and u = 0.6. 



4 Conclusion 



In this work we find that by applying a phase difference in the driving fields we 
can control chaos in Josephson junction. We also observe that for this type of 
control we have to choose the junction parameter values where the difference 
in voltages between the two junctions is negligible. We report that in a region 
where the difference in voltages between the two junction is not negligible 
chaos cannot be controlled by applying a phase difference in the driving fields. 
For the parameter region we choose the system was in periodic motion for 
a phase difference of ^ = 0.957r — l.Svr. Now we fix the phase difference as 
vr and vary other parameters such as dc bias, amplitude of applied field and 
coupling strength. It is found that even for large variation of these parameters, 
the system continues to be in periodic motion. So this might be of great 
practical importance as phase difference can be easily applied to the rf-field 
in an experimental set up. Thus it offers an easier way to control chaos and 
thus will provide an enhanced capability to design superconducting circuits in 
such a way as to maximize the advantages of non linearity while minimizing 
the possibility of instabilities. 
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Fig. 8. (a) shows the differential voltage plotted against time and (b) is the voltage of 
one junction with ^ = 0. It is observed that the variation is chaotic and the maximum 
difference in voltage is 0.05. In (c) and (d) shows the corresponding voltages with 
an applied phase difference of tt. Other parameter values are ag = 0.45, P = 0.15, 
iQ = 0.7 and u> = 0.6. 
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